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Abstract: 

In this paper the capacitated hub location problem is formulated by a minimax regret 
model, which takes into account uncertain setup cost and demand. We focus on hub 
location with multiple allocations as a strategic problem requiring one definite solution. 
Investigating how deterministic models may lead to sub-optimal solutions, we provide an 
efficient formulation method for the problem. 

A computational analysis is performed to investigate the impact of uncertainty on 
the location of hubs. The suggested model is also compared with an alternative method, 
seasonal optimization, in terms of efficiency and practicability. The results indicate the 
importance of incorporating stochasticity and variability of parameters in solving practical 
hub location problems. Applying our method to a case study derived from an industrial 
food production company, we solve a logistical problem involving seasonal demand and 
uncertainty. The solution yields a definite hub network configuration to be implemented 
throughout the planning horizon. 
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1 Introduction 


2 Literature Review 


Hub networks are one of the most common types of 
logistics systems serving urban transportation, airline 
networks, communication systems, and cargo networks. 
The basic characteristic of hub networks is that routing 
is performed through a subset of the connections between 
nodes instead of direct connections from the origin to 
the destination. The communication industry seems to 
be the first platform for using hubs, while decades has 
passed since hubs networks were institutionalized in 
logistics systems, transportation industry, air cargoes, 
and postal services. Nowadays, hub network design 
is a common practice for wholesalers, distribution 
companies, and food production industries whose main 
objective is to enhance logistics efficiency. Hub network 
configuration suggests using a set of hubs and spokes for 
connecting different origins and destinations. Different 
industries make use of hubs to deal with logistics 
activities in a more productive way by reducing direct 
transportation paths. Drawing an analogy might be 
helpful to clarify the importance of hub networks. 
A complete directed graph with K nodes has K x 
{K — 1) arcs, while all the nodes can be connected 
to each other by having a central node (hub) being 
connected to all the other peripheral nodes (spokes) 
which reduces the number of arcs to 2x(Ar —1). 
That is how the connectivity is achievable by utilizing 
fewer resources more productively. Hub location problem 
(HLP) originates from this idea; the challenge of deciding 
on allocation of hubs to obtain an efficient logistics 
network. 

The main objective in a common HLP is to 
minimize the total costs of establishing hubs and 
transportation of products between hubs and spokes. 
Hub location problems are categorized into capacitated 
and uncapacitated problems such that the former 
embodies most of the real-world problems. After HLP 
was introduced, subsequent problems like p-median, p- 
hub center, and hub covering problems were emerged to 
address different locational challenges of the industry. 

The principal purpose of p-median is to locate a 
number of hubs in the network so that the total 
transportation cost is minimized. The second problem, 
p-hub center, aims to optimize location of the hubs and 
allocation of nodes such that the cost of major routes in 
the network is minimized. In the hub covering problem, 
minimizing the total cost by finding the optimal location 
of the nodes and their corresponding allocation shapes 
the question where the number of hubs is not predefined. 
Such a problem introduces limits of coverage as the 
number of nodes that are connected to a hub is limited. 

Equally relevant to the problem type, main objective, 
and the decision variables are the questions of single 
and multiple network allocation patterns i.e. whether the 
spokes are to be connected to one hub or multiple hubs. 


HLP was first posed O’Kelly (19866|. The author 
introduced Single-HLP concerning assignment of the 
appropriate location to the hub and its connection to 
the spokes in a setting where there is no cost for hub 
establishment and it had infinite capacity. He then used 
a mathematical model to formalize another location 
problem on hub airline network O’Kelly (1987). The 
basic model made progress towards P-HLP, a quadratic 
model, incorporating a number of hubs with direct 
transportation routes Aykin ( 1990| . The proceeding 
linear model was extended to p-median location 
problem Campbell (1992) to incorporate multiple 
network allocation patterns. A more comprehensive 
mathematical model for multiple allocation HLP was 
presented Campbell (1994) embodying real world 
assumptions such as fixed cost for connecting spokes 
to hubs, minimum flow, and the capacity of nodes. 
Although most of the hub location models developed 
assumed potential locations for the hubs in a discrete 
space, there were early research studies relaxing this 
assumption and considering a continuous space |0’Kelly| 
(1986a), Aykin & Brown (1992). 

As the logistics operations got more complicated, 
new problems emerged with different objective functions, 
formulations, and solution techniques. The taxonomy of 
HLP suggests four ramifications of the original problems 
including Capacitated p-Median Problem [Shamsipour 


et al. 


(2012), HLP with star network structure Yaman 


(2008), Yaman & Elloumi (2012), p-hub center problem 
Campbell et al. (2005a|6), and p-hub covering problem 


Kara & Tansel (2003). For a more detailed review of 


the literature one may refer to Alumur & Kara (2008), 


Farahani et al. (2013). 


In the pioneering research studies the parameters 
were assumed to be deterministic though it was not 
realistic. Like the other location and logistics problems, 
incorporation of stochastic parameters is a promising 
research area that is receiving increasing attention. 
Some recent and more realistic approaches towards the 
model definition are as follows. pVIarianov fc Serra (2003) 
developed a model for locating hubs in the network of air 
transportation by formulating a M/D/c queuing system . 
The same research field was investigated by Yang (2009) 
who developed a two-stage stochastic programming 
model for air transportation with uncertainty in demand 
. They introduced a stochastic programming model to 
address the uncapacitated air freight hub location and 
flight routes planning under seasonal demand variations 


Yang (2009). 


Sim, Lowe, and Thomas, suggested a single 
assignment hub covering model where the arc travel 
time was normally distributed with given mean and 
standard deviation. The objective was to locate p hubs 
in order to minimize the longest transportation time 
in the network for a specified service level in delivery 
times. In their formulation, they did not include the 
cost of establishing the hubs and operating the network 
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I I Deterministic optimization 


I I Robust optimization 


Stochastic optimization 


I I Other approaches 


Figure 1 The evolution of optimization methods used for hub location problem: a transition from deterministic models to 
stochastic and robust models 


Sim et al. (2009). Moreover, Contreras et al. designed 


a model for multiple allocation HLP with uncertainty 
in both demands and transportation costs. They proved 
that stochastic problems with uncertain demands or 
dependent transportation costs are equivalent to a 
deterministic problem in which random variables are 
replaced by their expected values. However, in the 
case of uncertain independent transportation cost, the 
corresponding stochastic problem is not equivalent to 


the such a deterministic model Contreras et al. (2011). 


Snyder et al. presented a robustness measure that 
combines the two objectives by minimizing the expected 
cost while bounding the relative regret in each scenario 
for two classical facility location problems: the k-median 
problem and the uncapacitated fixed-charge location 


problem Snyder & Daskin (2006). In particular, their 


models seek a minimum expected cost solution that is 
p-robust. 


More recently, Alumur presented multiple allocation 
and single allocation HLP considering uncertainty 


factors in demand and setup costs Alumur et al.|(2012|. 


Abidi and Razmi suggested the same model as |Alumur| 
et al. (2012) (stochastic uncapacitated hub location 


problem with multiple allocations) for a similar case 
study where setup cost uncertainty was replaced by 
travel cost uncertainty [Adibi fc Razmi (2015). Some 
other approaches were suggested to move towards the 
reality in hub location network design. Risk pooling and 
its effects on hub location network design is suggested 


by Snyder et al. (2007). Reliability models for facility 


location are another approach taking the expected 
failure cost into account Snyder & Daskin (2005). 
Hult, Jiang, and Ralph solved the uncapacitated single 
allocation p-hub center problem with stochastic travel 
times using cutting planes and Benders decomposition 


Hult et al. (2014). For a complete review on facility 


location problems with different sources of uncertainty, 
Campbell & O’Kelly (2012). Figuredemonstrates 


see 


the evolution of methods used for solving hub location 
problems through time. A transition from deterministic 
optimization method to more sophisticated methods like 
robust and stochastic optimization is evident. 


As stated earlier hub location is an essential part of 
the strategic planning for distribution companies having 
far-reaching effects on their operational issues and 
productivity. Besides, as logistics activities are changing 
within time, the data used in hub network decisions 
can become outdated by the network utilization phase. 
Therefore, some of the parameters required for designing 
the network cannot be determined accurately. The most 
common uncertain parameters are costs, distances, and 
demands. Failure to consider uncertainty of parameters 
may lead to obtaining sub-optimal network designs as 
the determinant input parameters change. The cause 
of uncertainty of such parameters are as follows. The 
volatility of costs for initial procurement such as land, 
industrial equipment, and construction makes the setup 
cost uncertain. Although, the demand can be predicted 
by market research, the lag time between designing 
the network and its actual utilization makes many 
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Table 1 Basic model notation 

= {1, 2,n} Set of nodes 
Fk Fixed setup cost k G N 

dij Distance from node i G N to node j G N 

13 Collection cost 

5 Distribution cost 

a Transfer cost 

Wij Flow routed from node i G N to node j G N 

Cijkm The transportation cost from node i G N to node j G N routed via the hubs k G N and m G N 

Ffe The capacity of node k G N \i a, hub is located there 

Xijkm Fraction of flow from node i G N to node j G N routed via the hubs k G N and m G N 

yk Hub or spoke assignment for node k G N 


predictions outdated especially for the case of time- 
dependent demands like seasonal products. This issue 
indicates that uncertainty should also be considered for 
demand parameters. In many cases uncertain parameters 
follow a familiar probability distribution which needs 
stochastic optimization. There are also cases where 
the data does not fit to familiar distributions which 
require robust programming to take uncertainty into 
consideration. In both situations considering different 
scenarios in a discrete probabilistic space would take 


uncertainty into account Alumur et al. (2012). 


Our contribution: 

This research aims to investigate the effect of uncertainty 
on the solutions obtained from different modeling 
techniques proposed by the contemporary researchers. It 
concerns different approaches toward the formulation of 
a typical HLP with comprehensive uncertainty factors 
and capacity constraints. In many real world problems 
capacity constraints are an indispensable component of 
hub location mathematical modeling. 

However, to the best of our knowledge there was 
no published research studies on the investigation of 
uncertainty impacts on real hub location problems with 
capacity constraints and inaccurate decision parameters. 
This research study fulfills a needed gap on capacitated 
HLP under uncertainty as suggested for future research 
by many recent research studies Alumur et al. (2012), 


Adibi & Razmi (2015), Ghaffari-Nasab et al. (2015). In 


this study, a novel mathematical model is proposed to 
question this hypothesis whether deterministic modeling 
and seasonal optimization can be sound measures for 
obtaining the optimal location of the hubs under 
different sources of uncertainty. 

The structure of rest of the paper is as follows. Section 
presents a deterministic model to introduce the model 
foundation. Section suggests a more sophisticated 
optimization model with uncertain parameters. The 
efficiency of optimization model is evaluated in the 
Section by analyzing a numerical example. Finally, a 
practical case is discussed in Section to demonstrate 
the impact of the proposed approach in solving real world 
problems. 


3 Basic Model 

A review on the current literature reveals the presence 
of research studies on HLP with capacity constraints 


Ebery et al. 

( 

2000), 

Correia et al. 

( 

2011 

]2014 


However, the capacitated hub location literature requires 
realistic approaches towards incorporating uncertainty 
in modeling such as minimax regret model. This study 
suggests a novel modeling approach with four linear 
constraints to deal with uncertainty in demand and 
setup costs. As mentioned earlier, the foundation of 
optimization model is first introduced by replicating a 
deterministic model originally developed by | Campbell] 
(1994) and then the main approach of modeling is 


outlined in the next section to analyze capacitated HLP 
with multiple allocations. Adapting from a deterministic 
model proposed in a well-known study by | Campbell] 
(1994), the basic model notation is stated in Table]^ The 
total cost is structured in a way that it incorporates hubs 
setup cost as well as three types of transportation costs 
including collection cost for factory to hub, transfer cost 
for hub to hub, and distribution cost for hub-to-spoke. 
Accordingly, the transportation cost is formulated in Eq. 

0 - 


Oijkm — I3dik cydkm F ddn 


Vi, j, k,m G N (1) 


The capacitated HLP with multiple allocations is 
formulated as follows in ([^ to ([^. 

min EEEE E FkVk (2) 

i j k m k 

S.t. 


^ ^ ^ijkm E Vk 
m 

Vi, j,k G N 

(3) 

^ ^ ^ijkm — Um 
k 

1 Vi, j, m G N 

(4) 

i j m 

^ijkm — ^kUk Wk G iN 

(5) 

Vk € {0,1} 

\/k G N 

(6) 

0 E ^ijkm E 1 

Vi, j, k,m G N 

(7) 
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Table 2 Minimax regret model notation 

S'j: All the scenarios with different uncertain setup costs s' € S'j: 

F^, The cost of establishing a hub at node k in scenario s' 

Sw All the scenarios with different uncertain demands s € 

Ps The probability that scenario s G occurs 

Wfj Demand routed from node i to node j in scenario s G Sw 


In this formulation, Xijkm stands for the decision variable 
which is the fraction of flow from location i (origin) 
to location j (destination) that is routed via hubs at 
locations k and m in that order, yk represents a binary 
variable showing whether node fc is a hub by taking one, 
or it is a spoke by taking zero. Eq. (§ and a ensures 
that the flow passes through the hubs. Eq. ([^ takes 
the capacity of the hubs into account. Einally, domain 
constraints are formulated in ^ and 0. 


4 Minimax Regret Model 


The impact of hub location on total cost makes it a 
crucial strategic decision to be made in logistics. As an 
attempt to deal with lack of precise information on the 
operational parameters of logistics networks, minimax 
regret model can be deployed. Minimax regret is a 
robust optimization approach to minimize the worst-case 
regret. The aim of this technique is to perform as closely 
as possible to the optimal course. Since the minimax 
criterion applied here is to the regret rather than to 
the payoff itself, it is not as pessimistic as the common 
minimax approach. 

One benefit of minimax regret models are that they 
are independent to the probabilities of the various 
outcomes. Morover, in comparison with other methods 
of robust optimization, it is easy to use specially 
in cases where regret can be accurately computed 
while probabilities of scenarios are difhcult to estimate. 
Although minimization of maximum regret is a common 
practice, many other techniques have been proposed 
for both stochastic and robust optimization in HLP. 
Ghaffari-Nasab et al. considered the capacitated single 
and multiple allocation HLP with stochastic demands 


Ghaffari-Nasab et al. (20151. The main emphasis 


of their paper relies on robust optimization tools 


developed by Bertsimas & Sim (20041 for linear 


programming problems. Situation and strategies in 
logistics decision making need different approaches to 
tackle the uncertainty. Eor a review of the application of 
such approaches to facility location problems, see | Owen] 
& Daskinl (11998 ) and ISnyderl (p006|) . 


It is assumed in this paper that the uncertainty 
in demand can be described by considering a limited 
number of scenarios. It is also assumed that in 
each scenario demand parameters are certain values. 
Moreover, uncertain behavior of setup costs is assumed 
to be interpretable by considering different scenarios. 
Deploying such scenarios alongside minimax regret 


programming, the model will be able to tackle real world 
problems in an uncertain environment Alumur et al.| 


(2012). The minimax regret model notation as originally 


suggested by Alumur et al. (2012) for another type of 
HLP is stated in Table |2j 


Z*, = min WfjCijkmXijkm 

s i j k m 


( 8 ) 


S.t. 


^ ^ ^ijkm ^ Vk 

771 

Vz, j,k G N 

(9) 

^ ^ ^ijkm ^ ym 
k 

1 Vz, j, m G N 

(10) 


Xijkm ^ ^ kUk V/z, S G Al 

(11) 

i j 771 



Vk € {0,1} 

\/k G N 

(12) 

0 ^ijkrTi ^ 1 

Vz,_), k,m G N 

(13) 


Z*, is the optimal solution of the model above. 
Respectively as the exact scenario that will occur is 
not known, a minmax regret model can be considered 
as in Eq. to Eq. in which the maximum regret 
is to be minimized. In this mathematical model, the 
representation of yk is the same as before. The above 
model can be easily linearized by defining variable R 


such that R> R'g V s' G S'j. 


min max = R' 

(14) 

S.t. 


Xijkm < Vk Vz, j,kG N 

771 

(15) 

Xijkm < Um Vz, j,m G N 

k 

(16) 

< r kVk yk,s gn 

(17) 


i j m 


Rs — Ps WijC-ijkmXijkm + 

S i j k 771 

J2F^,yk-Z:, Ws'GS'f 

k 

yk G {0,1} WkGN 

0 ^ Xijkm ^ 


(18) 

(19) 

( 20 ) 
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Table 3 Test case 1: different solutions for capacitated HLP with multiple allocations 



a=0.3 


a=0.5 


P 

II 

o 


a=l 



Gost 

Hub 

Cost 

Hub 

Cost 

Hub 

Cost 

Hub 

BDM 

2905117 

2,3 

2989450 

1,3 

3065952 

1,3 

3138530 

1,3 

s/i 

2884970 

2,3 

2969830 

2,3 

3054288 

2,3 

3138440 

2,3 

s/2 

2084747 

2,3 

2169607 

2,3 

2254065 

2,3 

2338217 

2,3 

«/3 

2461068 

2,4 

2547230 

2,4 

2630818 

2,4 

2712427 

2,4 

S/4 

1779440 

1,3 

1864380 

1,3 

1942708 

1,3 

2018130 

1,3 

MRM 

- 

3,4 

- 

3,4 

- 

3,4 

- 

3,4 


5 Computational Analysis 

To test the model proposed in Section some test 
problems associated with the distribution of products 
among five different cities 0-lE were used (highlighted 
cells to be discussed later). The other parameters were 
considered as (5 = 5=1 and a € {0.3,0.5, 0.7,1}. Four 
different scenarios for uncertainty in the setup cost were 
designed in which F^, was randomly selected. Moreover, 
for analyzing the uncertainty in demands, four different 
scenarios were selected with equal probabilities Pg = 
0.25. GAMS software was used to solve the numerical 
example. 

Optimal solutions of different modes are represented 
in Table m The basic deterministic model is abbreviated 
to BDM and is shown in the first row. The four scenarios 
of the stochastic model are presented in the next rows 
according to four different values assumed for cost of 
land represented by s/. Finally, minimax regret model 
is abbreviated to MRM and is demonstrated in the 
lowest row. The problems are solved separately based 
on each scenario considering minmax regret models. The 
solutions can be compared with that obtained from 
the basic deterministic model in which setup costs and 
demands were set to the mean values. 

The second column form left side of Table|3]shows the 
total cost and the third column represents the optimal 
location of the hubs. Note that as the objective function 
in the minmax regret model is represented differently, 
total cost cannot be compared with that of scenarios. 
Hence, the cells for values of the minimax regret cost 
are left empty. As was evident in Table the optimum 
location of the hubs in the minimax regret model differs 
from the other scenarios. It can be concluded from the 
observation that it is more appropriate to use a minimax 
regret model instead of estimating costs and demands or 
using a deterministic scenario. 

No relationship is observed between the costs of 
setting up a hub and selecting a location for the hub. 
For example, node 4 has the highest setup cost, but 
in some problems it is selected as the optimal location 
of hub. This observation indicates that in addition to 
the setup cost, the demand and geographical location 
are also determinant factors in the optimal solution of 
capacitated HLP. 

In order to show the robustness of the model to re¬ 
designing the position of the points, a new city is added 
to the test case. The highlighted rows and columns in 


Tables[^-[^show the added city distances and cost under 
different scenarios as well as demand from the other 
cities. The modified problem is solved using GAMS. 
As expected, the optimal location of hubs for MRM 
would not change from cities 3 and 4 which shows the 
robustness of minimax regret model to points of sale. 
Furthermore, the minimax regret model is robust to the 
variation of demand and cost as already demonstrated 
in Table [S] 

As an attempt to investigate larger test cases, 
randomly generated problems with 5,6,10 and 15 nodes 
are analyzed in which (i = 5 = 1 and a = 0.5. For each 
problem, 4 scenarios of setup cost and 4 scenarios 
of demand (with equal probabilities Pg = 0.25) are 
considered. The results are reported in Table As 
it is illustrated in Table the optimal location of 
hubs can be changed based on the modeling method 
(deterministic versus minimax regret) and variations 
in demand. This supports the previous argument that 
change of parameters requires re-configuration of hub 
network if the model is deterministic. However, taking 
the variations of parameters into account, a minimax 
regret model can be used to determine a definite hub 
network configuration that works well under uncertainty. 
This approach is implemented on a practical case in the 
next Section. 

Table 4 The optimal location of hubs for the other test 
cases 


Test case 

2 

3 

4 

5 

No. of nodes 

5 

6 

10 

15 

BDM 

1,2 

1,3 

5,3 

5,2 

s/i 

3,4 

1,3 

3,9 

5,11 

s /2 

1,4 

1,3 

1,5 

2,4 

s/3 

1,3 

1,4 

2,3 

1,5 

s/4 

1,3 

1,4 

1,9 

3,13 

MRM 

1,3 

1,3 

5,9 

1,12 
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6 Practical Case Study 


Table 5 Demands in different seasons 


The application of proposed model is evaluated using 
the data from an Iranian industrial food production 
company. The case of chocolate production in Shirin Asal 
Tabriz Co. is a well-known local hub location example, 
investigated by Rostami et al. (20121. They analyzed 
the impact of estimated demand on the configuration 
of hub networks in different scenarios. According to 
their model, the location of hubs can be changed 
seasonally. In contrast, they did not consider the costs 
of seasonal re-configuration the network as the model 
required so. As already discussed, we consider HLP 
as a strategic decision making process requiring an 
unchangeable solution as proposed by this research. 

As mentioned earlier, the fluctuation of prices makes 
hub establishment an uncertain activity in terms of 
monetary issues. Moreover, according to the data 
reported by the company, the demand for products is 
seasonal and the setup cost is highly dependent on time, 
making it an appropriate case to be analyzed by our 
model due to the fluctuations in demand and variability 
of the setup costs. The scenarios are designed by dividing 
the year into four seasons with equal length for demand 
and considering five scenarios for setup costs with O.TF^. 
to l.SPfc. 

The geographical structure of Shirin Asal Tabriz Co. 
market is as follows. The main factory is located in 
the city of Tabriz in the north-west of Iran supporting 
36 distribution points with demand across the country. 
The national market can be divided into three regions 
namely: the west, the center and the east. The data 
presented in this study are related to the demand of the 
west part of the country with 14 nodes. The company 
management sought to establish hubs among these cities 
(14 locations). Table outlines the demands in each 
scenario. Table includes the capacities and the cost (in 
million Rials) of establishing hubs in each location. The 
problem is to find the best assignment of the hubs to the 
cities and their allocation according to nodes capacities 
and uncertain demand and setup cost. According to 
the studies performed for this particular example, a = 
0.4,/3 = 1 and ^ = 1 are calculated [Rostami et al.| 
(20121 and Pg = 0-25 to represent four seasons. The 
optimization is performed by GAMS software on a 
personal computer with 2.4 GHz Intel Core i5 CPU and 
4 GB of RAM. 


To demonstrate the practical application of the 
proposed minimax regret model, a deterministic model 
(setup costs and demands set to the mean values) is 
first solved and then the minimax regret model is used 
to analyze the case with respect to significant sources 
of uncertainties discussed. The result of deterministic 
model and minimax regret model are shown in Fig. 
where the flow between the factory and the hubs (cost 
represented by /3) as well as that between the hubs (cost 
represented by a) are represented by thick lines and the 
flow between hubs and spokes (cost represented by <5) 
is represented by narrow lines. The network obtained 


City 

Spring 

Summer 

Fall 

Winter 

Rasht 

9205 

7899 

21848 

26510 

Kermanshah 

10459 

11256 

18751 

19628 

Tabriz 

48022 

39529 

65890 

94831 

Tehran 

14412 

12571 

44070 

91906 

Zanjan 

10590 

11732 

20402 

21218 

Qazvin 

4424 

5995 

10848 

12502 

Hamedan 

6270 

4802 

9273 

9505 

Urmia 

17022 

16006 

24951 

23234 

Ardabil 

13764 

19839 

19281 

16767 

Sanandaj 

7996 

6105 

11330 

10807 

Shahrekord 

6142 

5721 

10320 

9065 

Ham 

4044 

4135 

6689 

7881 

Karaj 

18519 

22050 

41018 

48957 

Arak 

4272 

3726 

9402 

10287 


Table 6 Capacity and setup cost for different cities 

City 

Capacity 


Setup cost 


Min 

Avg. 

Max 

Rasht 

275200 

910.56 

1300.80 

i691.04 

Kermanshah 

154200 

431.76 

616.80 

801.84 

Tabriz 

321640 

900.59 

1286.56 

1672.52 

Tehran 

38528 

1358.78 

1941.12 

2523.45 

Zanjan 

201360 

703.80 

1005.44 

1307.07 

Qazvin 

123840 

346.75 

495.36 

643.96 

Hamedan 

97600 

329.28 

470.40 

611.52 

Urmia 

110080 

308.22 

440.32 

572.41 

Ardabil 

123840 

346.75 

495.36 

643.96 

Sanandaj 

103200 

288.96 

412.80 

536.64 

Shahrekord 

137600 

385.28 

550.40 

715.52 

Ham 

82560 

231.16 

330.24 

429.31 

Karaj 

460960 

1500.68 

2143.84 

2786.99 

Arak 

82835 

231.93 

331.34 

430.74 


from the deterministic model is shown on the map on 
top of Fig. [^suggesting establishing two hubs in cities of 
Ardebil and Kermanshah, basing the decision on mean 
values and ignoring the uncertainty. On the other hand, 
the minimax regret model suggests establishing hubs in 
three cities of Qazvin, Zanjan and Arak as the solution 
is illustrated on the map in the bottom of Fig. 

It is noteworthy that in the optimal solution obtained 
by the minimax regret model, two hubs located in cities 
of Qazvin and Zanjan are connected. Such solutions 
are not likely be optimal in an uncapacitated setting. 
However, in capacitated problems if high demand nodes 
are close to low capacity nodes with low setup cost, we 
may observe the optimal solution connecting hubs. In 
this case, the connectivity of hubs in Qazvin and Zanjan 
is mainly due to the high demand in cities of Tehran and 
Karaj and relatively low capacity of Qazvin as their close 
hub with low setup cost. 

Furthermore, the most striking observations to 
emerge from the network configuration are the 
distribution routes from the company to cities of Ardebil 
and Urmia. As evident in Fig.j^the products need to be 
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first collected in Zanjan (hub) and then be distributed 
to the two cities of Ardebil and Urmia. Although based 
on the distances such an optimal solution might seem 
odd, the hub establishment cost values render this 
reasonable. Direct transportation of products from the 
factory requires Tabriz to be a hub, while Tabriz is the 
fourth most expensive locations for hub establishment. 
Worthy of mention is that the factory cannot support the 
functions of the hub and for the purposes of the problem 
Tabriz is similar to the other cities. 



Figure 2 Different mathematically optimal hub location 
solutions obtained from deterministic model (top) 
minimax regret model (bottom) 



<P 

Figure 3 Comparison of total cost from seasonal 

optimization and minimax regret model based on 
different values of (j> 


time Rostami et al. (2012). The reported cost is exclusive 
of network seasonal re-configuration cost involved in 
operations dealing with training, human resources, and 
facility transportation. If we consider (j) as the ratio 
of re-configuration cost (compromising of operations 
dealing with training, human resources, and facility 
transportation) to the total cost of hub network 
establishment, results from the two models can be 
compared as in Fig. showing a break-even point for 1 
year operation. As evident in the figure, our suggested 
model is capable of reducing the total cost for </> > 
0.0174 and the cost reduction is proportional to the 
ratio defined 4>- Therefore, for industries where hub 
network is easy to be re-configured every season based 
on demand variations, seasonal optimization suggested 
by [Rostami et ah (2012) can be used. For the other 
cases where the company needs to spend a non-trivial 
{(j) > 0.0174 in here) amount of money on hub network 
re-configurations, a definite hub network design as 
suggested in the present study is more cost-effective. 

The substantial difference between the solutions 
obtained from the deterministic model and minimax 
regret model shows the impact of uncertainty on 
HLP. The differences in hub network designs let us 
argue that deterministic analysis of HLP suffers from 
capability of industrial practice in cases where there 
are source of uncertainty. Deterministic analysis may 
lead to sub-optimal solutions imposing excessive costs to 
the company in the long term. Therefore, deploying a 
definite solution obtained from a comprehensive model 
is suggested to deal with logistics with seasonal demand 
and uncertainty to be implemented throughout the 
planning horizon. 


The total cost of the network obtained by seasonal 
optimization model is computed a 960 day operation 
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7 Application Perspectives 

Designing logistics network under uncertainty has 
applications already in demand in various businesses 
across different industries including those with a seasonal 
demand and other time-variant decision parameters. 
This study aims to broaden the hub location body 
of knowledge via a quantitative investigation of an 
industrial problem to assist executives by providing 
a novel network design model. The suggested model 
is to be used for increasing cost-effectiveness as 
a competitive advantage as well as meeting other 
operational objectives. 

In more specific terms, this study suggests a method 
to tackle uncertainty in locational decision making 
to meet essential technical challenges associated with 
designing hub networks like inaccuracy of estimates and 
variability of demand. A minimax regret model can be 
used in decision making problems where the uncertainty 
can be modeled by considering a number of scenarios 
for the uncertain parameter while no historical data is 
available to estimate probability distributions required 
for stochastic programming. Practical examples provided 
show how the suggested technique can be used to obtain 
a definite logistics network under uncertain demand and 
setup cost. 

8 Conclusion 

In this paper, multiple allocation capacitated HLP was 
investigated in a setting where setup cost and demand 
were both uncertain. After outlining a well-know model 
in Section different extensions are discussed and a 
minimax regret model was developed in Section for 
considering such sources of uncertainty. A computational 
analysis was performed in Section to investigate the 
changes in the optimal location of the hub caused 
by different modeling approaches. The result showed 
that the optimal solution changes where the model is 
associated with uncertain parameters. According to the 
numerical example in Section]^ ignoring the uncertainty 
may change the whole hub location solution drastically. 
Moreover, the industrial application of the proposed 
method was addressed to by discussing a case from 
an industrial food production company. The results 
obtained from solving the minimax regret model for the 
case study and comparing it to the alternative modeling 
technique confirmed the efficiency of our suggested 
model in providing a definite solution for an industrial¬ 
sized problem challenged by an uncertain environment. 

There are a number of research directions that are 
hoped to be investigated in the near future. Firstly, 
one may introduce a more sophisticated problem with 
pervasive sources of uncertainty by parameters not 
necessarily following a familiar distribution function to 
be solved by stochastic programming and compared with 
the technique used in this paper. It would contribute 
to our idea if the research is associated with a critical 


industrial-sized case of an uncertain environment and 
counterproductive practices. Secondly, it is suggested for 
further research to tackle large HLPs of this type by 
tailored evolutionary algorithms and analyze practical 
cases for different datasets such as international couriers. 


A Computational analysis and case study 
data 

Table 7 Distance between 5 cities used in the 
computational analysis example 


City 

1 

2 

3 

4 

5 

1 

0 

590 

485 

325 

348 

2 


0 

588 

526 

414 

3 



0 

599 

280 

4 




0 

319 

5 





0 

6 

729 

456 

578 

423 

253 


Table 8 Four scenarios of demand used in the 
computational analysis example 


Demand scenario 1 

City 

1 

2 

3 

4 

5 

1 


12612 

18085 

16652 

63467 

2 

86249 


92152 

42150 

6979 

3 

15295 

53550 


87155 

81086 

4 

86940 

90960 

47946 


88818 

5 

61337 

91633 

76636 

91141 


6 

45400 

54222 

35200 

48576 

43253 

Demand scenario 2 

City 

1 

2 

3 

4 

5 

1 


68051 

78747 

43698 

48346 

2 

71429 


67029 

38488 

44321 

3 

39460 

28955 


73469 

62608 

4 

63443 

7932 

90296 


68353 

5 

19322 

12575 

6864 

20751 


6 

54378 

68782 

42200 

78245 

65200 

Demand scenario 3 

City 

1 

2 

3 

4 

5 

1 


49129 

72170 

91123 

80320 

2 

65650 


26966 

53580 

26892 

3 

63409 

34737 


16355 

77912 

4 

18540 

57047 

67416 


25912 

5 

14567 

24116 

84892 

27186 


6 

86245 

9237 

83240 

9257 

8424 

Demand scenario 4 

City 

1 

2 

3 

4 

5 

1 


35764 

29767 

10636 

15561 

2 

21637 


72711 

8641 

55549 

3 

26601 

57047 


52084 

46490 

4 

59851 

53809 

38386 


4810 

5 

46845 

87287 

55457 

88819 


6 

16784 

17865 

12232 

18526 

15423 
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Table 9 Capacities and four scenarios of setup cost used in the computational analysis example 


City 

1 

2 

3 

4 

5 

6 

Capacity 

682423 

765892 

876543 

986578 

546879 

454786 

Setup cost in scenario 1 

1414016725 

965058131 

1483010032 

1529540551 

1607592118 

1825643271 

Setup cost in scenario 2 

1213461250 

727806327 

920038779 

1826350935 

818597258 

1934765489 

Setup cost in scenario 3 

1710445940 

1329737654 

1936001128 

720329416 

1203047686 

1745980390 

Setup cost in scenario 4 

758042396 

1890660867 

622937195 

1643029427 

1699171001 

1693049714 

Average setup cost 

1273991578 

1228315745 

1240496784 

1429812582 

1332102016 

1799859716 


Table 10 The distance between 14 cities of the practical case study 



1 2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

1. Rasht 

0 590 

485 

325 

348 

185 

401 

739 

266 

565 

868 

774 

285 

577 

2. Kermanshah 

0 

588 

526 

414 

433 

189 

582 

791 

136 

731 

184 

538 

365 

3. Tabriz 


0 

599 

280 

455 

609 

308 

219 

452 

1142 

772 

574 

785 

4. Tehran 



0 

319 

150 

337 

907 

591 

501 

543 

710 

50 

239 

5. Zanjan 




0 

175 

329 

588 

377 

278 

862 

598 

282 

505 

6. Qazvin 





0 

244 

763 

451 

453 

584 

617 

106 

303 

7. Hamedan 






0 

610 

667 

164 

568 

373 

354 

176 

8. Urmia 







0 

527 

446 

1178 

766 

729 

786 

9. Ardabil 








0 

655 

1134 

975 

552 

843 

10. Sanandaj 









0 

732 

320 

523 

340 

11. Shahrekord 










0 

719 

579 

392 

12. Ham 











0 

706 

514 

13. Kara) 












0 

322 

14. Arak 













0 
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